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AN INVESTIGATION OF B LEED CONFIGURATIONS AND THEIR EFFECT 
ON SHOCK WA VE/BOUNDARY LAYER INTERACTIONS 


ABSTRACT 

The design of high efficiency supersonic inlets is a complex task involving the 
optimization of a number of performance parameters such as pressure recovery, spillage, 
drag, and exit distortion profile, over the flight Mach number range. Computational 
techniques must be capable of accurately simulating the physics of shock/boundaiy layer 
interactions, secondary comer flows, flow separation, and bleed if they are to be useful in 
the design. In particular, bleed and flow separation, play an important role in inlet unstart, 
and the associated pressure oscillations. Numerical simulations were conducted to 
investigate some of the basic physical phenomena associated with bleed in oblique shock 
wave boundary layer interactions that affect the inlet performance. 


INTRODUCTION 

The experimental measurements in most mixed compression inlet tests are limited to 
the static pressures over the inner surfaces and total pressure recovery at the e ngin e face 
[1]. In a few experimental studies, some boundary layer profile measurements were 
obtained from pitot pressure rakes ahead and behind the shock/boundary layer/bleed regions 
[2-5], Comparisons of flow computational results with the experimental measurements in 
supersonic inlets revealed discrepancies in the prediction of shock locations, and velocity 
profiles after the interactions [6-8], A better understanding of these interactions is essential 
to improving inlet performance predictions.. 

The purpose of the present investigation is to characterize the flow field in the 
oblique shock wave/boundary layer/bleed interaction region and to develop a basic 
understanding of flow separation control mechanisms. The flow computations were carried 
out for different bleed configurations, and shock strengths. The solution domain extended 
across the interaction, and inside the bleed slots. A systematic investigation was conducted 
to determine the effects of the bleed location relative to the shock, slot width, depth and slot 
angle, on the flow characteristics and the bleed performance parameters at different bleed 
mass flow rates [9-13]. The obtained flow characteristics inside the bleed port were 
Compared with the recently obtained experimental measurements [14, 15] with and without 
incident shock at NASA Lewis Research Center. 

The results of the bleed investigation were presented in two formats, one giving the 
flow details in order to understand the nature of the interaction, and another giving the 
bleed performance parameters such as the discharge coefficient and total pressure recovery 
which are useful to the designer. The detailed pressure and Mach number contours as well 
as velocity profiles within the interaction zone and inside the bleed slots revealed the 


important flow physics and the mechanisms affecting the flow separation and demonstrated 
the variation in the bleed flow velocities, and pressure gradient distribution across the bleed 
port. On the other hand, the bleed performance results in the form of bleed efficiency and 
pressure recovery at different plenum pressures and for different bleed configurations is 
useful in selecting the optimum bleed configuration in inlet designs. 


ANALYSIS 

The solution domain for the compressible Navier-Stokes equations extended inside 
the bleed slot and upstream and downstream the interaction region. The change in bleed 
(suction) mass flow was accomplished in the numerical investigation throu gh changing the 
static pressure at the bottom of the bleed slot over a wide range covering choked and 
unchoked flow conditions. The PARC code [16] was selected after it was validated [10] 
for the case of separated SWTBLI at Mach 2.96 through comparison with Law’ .5 [17] 
experimental data. The mesh spacing was varied in both the streamwise, and normal 
directions to capture the important flow characteristics in the shock wave boundary layer 
interactions, and inside the bleed port. 


SUMMARY 


I- Eff^t of Normal Bleed Slot Location in Lami nar Boundary Laver/Shock 

Interaction 

A study was conducted for bleed through a normal slot at three locations upstream, 
across and downstream the shock [9]. The external flow conditions were M® — 2.0, 

Re® = 1.812xl0 6 /ft., and an impingement oblique shock angle of 32.585 deg., 
corresponding to the experimental test conditions of Hakinen et al. [18] without bleed. The 
bleed slot width of D = 0.01 1 inch was equal to the boundary layer displacement thickness 
5*before the interaction at x = 1 .96 inch from the flat plate leading edge (D/5*=l) and the 
slot depth was three times its width (L/D = 3.0). For a given plenum pressure at the slot 
exit, it was found that the shock induced flow separation on the plate surface was eliminated 
in the across the shock bleed configuration, but only reduced for the other two bleed 
locations. The results of this investigation were published in the Journal of Propulsion and 
Power, Vol. 11, No. 1, pp. 42-48, Jan. -Feb. 1995, and presented as AIAA Paper 91-2014 
entitled “The Effect of Bleed Configuration on Shock Laminar Boundary Layer 
Interaction,” at the Joint Propulsion Conference in Pasadena, CA, June 14-16, 1991 
(Appendix A). 


2 . 


Effect of Bleed Mass Flow Through 
Boundarv-Laver Interaction 


Wl'lMU 


al Slot in Shock-Wave/Turbulent- 


A study was completed to investigate the effect of bleed mass flow in a normal slot 
across the incident oblique shock impingement point. The external flow conditions were 


M® - 2.96, Re® = 1.2xl0 7 /ft., and impinging oblique shock angle of 25.84 degree, 
corresponding to the experimental test conditions of Law [17] without bleed. The bleed 
slot width was D = 0.0213 ft. which is equal to 1.617 times the boundary layer thickness 6 
upstream of the interaction at x = 0.9 ft. (D/8 = 1.617), and the bleed mass flow varied up 
to choked conditions (30% of the incoming boundary layer mass flow). The computed 
bleed results [10] indicated a complex expansion and compression wave structure across the 
slot opening due to bleed flow entrainment. The net effect of the bleed on the flow 
downstream is to reduce the momentum and displacement thickness and to move the 
reflected shock closer to the surface. The reduction in the momentum thickness 
downstream of the interaction increased with the increase in bleed mass flow, but eventually 
plateaud. This work was published in the Journal of Propulsion and Power. Vol. 10, No. 1, 
pp. 79-81, Jan.-Feb. 1995, and presented as AIAA Paper 92-3085 entitled “An 
Investigation of Shock/Turbulent Boundary Layer/Bleed Interaction,” at the Joint 
Propulsion Conference in Nashville, Tennessee, July 6-8, 1992 (Appendix B). 

3. Effect of Normal Bleed Configuration 

A parametric study was conducted to determine the effect of normal bleed slot 
configuration on shock wave turbulent boundary layer interactions [11]. The external flow 
conditions were M« = 2.96, Re® =1.2x1 0 7 /ft. and oblique shock impingement angle of 
25.84 degree. Different slot width, depth and inclination were investigated. The results of 
the parametric investigation of the normal bleed slot configuration in STBLI were presented 
as AIAA Paper 93-0294 entitled “A Parametric Study of Bleed in Shock Wave Boundary 
Layer Interactions,” at the 3 1st Aerospace Sciences Meeting in Reno, Nevada, January 1 1- 
14, 1993 (Appendix C). 

(a) Slot Width 

Flow computational results were obtained over a wide range of bleed mass flows for 
three different slot widths 0.8, 1.6 and 3.2 times the boundary layer thickness upstream of. 
the interaction [14]. The results indicate that the bleed discharge coefficient and total 
pressure recovery are not very sensitive to the bleed slot width. 

(b) Slot Aspect Ratio 

Two slot aspect ratios = 3 and 6 were investigated for the intermediate slot width. 
Both the discharge coefficient and the momentum thickness downstream of the interaction 
were sensitive to changes in L/D and both improved as the aspect ratio increased.. 

(c) gMfoqrtfol 

Additional flow computations were performed for three different bleed locations 
relative to the shock impingement point for the widest bleed slot (D/8 = 3.2) over a range of 
bleed mass flows up to choked conditions. The upstream, across, and downstream bleed 
slot locations corresponded to the slot center at -D, 0 and +D relative to the shock 
impingement point. The results indicate that across the shock bleed gives the largest 
reduction in the momentum thickness downstream of the interaction at all bleed mass flows. 


However, the total pressure recovery at the bleed slot exit was higher for downstream 
bleed. 

4. Effect of Bleed Slot Angle and Location 

Detailed flow computations [12] were performed for 20° slanted and 90° normal 
slot3 at three different locations upstream, across and downstream of the inviscid shock 
impingement location relative to the flat plate surface at Ho = 2.96, Re» = 1.2x107ft and 
oblique shock impingement angle of 25.84 deg. The slot opening was kept the same for 
both slots with D/5 = 0.8085 for the normal slot and 0.2766 for the 20° slanted slot. The 
results indicate a large recirculating flow region inside the normal slot but not in the 20° 
slanted slot causing the choked discharge coefficient to more than triple. Bleed through the 
slanted slot produced redeveloping boundary layer downstream, with a higher friction 
coefficient. The largest reduction in the boundary layer momentum and displacement 
thickness do wnstream of the interaction were associated with the normal slot bleed across 
the shock. 

The slanted bleed discharge coefficient was found to be more sensitive to bleed 
location and had the highest value when the bleed slot was located upstream of the shock. 
Normal bleed across the shock resulted in the largest reductions in both boundary layer 
displacement and momentum thickness downstream of the interaction. The results of this 
investigation will appear in Journal of Propulsion and Power, Vol. 1 1, No. 6, Nov.-Dee.. 
1995 entitled “Shock-Wave Boundary-Layer Interactions with Bleed, Part 2: Effect of Slot 
Location,” and were presented as AIAA Paper 93-2992 entitled “An Investigation of Shock 
Wave Turbulent Boundary Layer Interaction with Bleed Through Slanted Slots,” in 
Orlando, Florida, July 1993 (Appendix D). 

5. Plenum Interactions 

The numerical investigation of shock-wave/turbulent boundary-layer/bleed 
interactions was conducted with a plenum at the end of the bleed slot [13]. The results 
were compared with the experimental data obtained by Willis et al. [14] and Davis et al. 

[15] in the l*xl* supersonic tunnel at NASA Lewis Research Center. In this study the 
external flow conditions were M» = 2.46, Re* = 1.81xl0 7 /ft., and the ratio of the slot 
width was one third the incoming boundary layer thickness before the interaction (D/5 = 
1/3). The flow simulations were performed over a range of plenum pressures to simu late 
increased bleed mass flow rates up to choking for two shock strengths corresponding to 
shock generator wedge angle of 0° and 8°. The computational results with and without the 
plenum indicated that the interaction between the plenum and the recirc ulatin g flow inside 
the slot caused fluctuations in the bleed mass flow at the slot exit. These variations were 
associated with the changes in the size of the separated flow regions over both slot walls at 
exit. However, the bleed mass flow at the slot opening did not fluctuate and r emain^ 
steady at the mean value of the exit mass flow. The computed bleed mass flow at the slot 
inlet were in good agreement with the experimental measurements. These results were 
presented in AIAA Pacer 95-0033. "Flow Characteristics is in Boundary Layer Bleed Slots 


4 


with Plenum,” at the 33rd Aerospace Sciences Meeting, in Reno, Nevada, January 9-12, 
1995 (Appendix E). 
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Effect of Bleed Configuration on Shock/Laminar 
Boundary-Layer Interactions 

A. Harried* 

University of Cincinnati, Cincinnati, Ohio 45221 
and 

T. Lehnigt 

Cold Jet, Inc., Loveland, Ohio 45140 


The fiowfield characteristics are simulated numerically in an oblique shock wave/laminar boundary-layer 
interaction for three different bleed configurations. The numerical solution for the fiowfield is obtained for the 
strong cOnservation-law form of the two-dimensional compressible Navier-Stokes equations using an implicit 
scheme. The computations model the flow in the Interaction region and inside the bleed slot tor an impinging 
oblique shock of sufficient strength to cause boundary-layer separation on a flat plate boundary layer. The 
computed results are presented for a normal bleed slot at three different locations: I) upstream. 2t across, and 
3) downstream of the shock impingement point. The detailed flow characteristics in the interaction zone and 
inside the bleed slot are compared for the three cases. The resulting surface pressure and shear stress and 
boundary -layer displacement and momentum thickness distributions are also compared for the three Weed cases 
and for the shock wave/boundary-layer interaction without bleed. 


Introduction 

A IR bleed systems are used for controlling the shock wave/ 
boundary-layer interactions when operating at super- 
sonic speeds. Proper bleed system design is particularly im- 
portant for the efficient and stable operation of mixed 
compression supersonic inlets. 1 The fundamental objectives 
of supersonic inlet bleed system design are to provide good 
aerodynamic flow characteristics with minimum boundary- 
layer bleed. 5 -* The elimination or reduction of separation in 
shock wave/boundary-layer interactions 4 is essential for con- 
trolling the pressure loss and maintaining a stable fiowfield 
in supersonic inlets. Minimizing the bleed required to achieve 
boundary-layer control allows the penalties associated with 
incorporation of a supersonic inlet bleed System to be mini- 
mized as well. These penalties include a loss of inlet mass 
flow and an increase in drag due to venting of the bleed air 
into the freestream. Ideally, removal of only that portion of 
the low-momentum boundary layer needed to inhibit sepa- 
ration is desired. 

Comparisons of internal flow computational results' " 7 with 
the experimental measurements in supersonic inlets 4 " re- 
vealed reasonable agreement between the computed and 
measured surface pressures upstream of the ramp bleed. How- 
ever. discrepancies in the predicted shock locations and ve- 
locity profiles were observed downstream of shock/boundary- 
layer interactions with bleed. 

Hamed and Shang 1 reviewed the existing experimental data 
base for shock wave/boundary-layer interactions in supersonic 
inlets and other related configurations. According to this sur- 
vey. there is enough experimental evidence 1 '" 14 to indicate 
that local bleed can control flow separation in shock wave/ 
boundary-layer interactions. There are disagreements, 1 how- 
ever, among the different experimental results regarding the 

Presented as Paper 91-2014 at the AlAA/SAE/ASME/ASEE 27th 
Joint Propulsion Conference, Sacramento, CA. June 24-26. 1991; 
received Aug. 2. 1991; revision received May 4. 1993; accepted for 
publication May 9. 1994. Copyright © 1991 by the American Institute 
of Aeronautics and Astronautics, Inc. All rights reserved. 

•Professor. Department of Aerospace Engineering and Engineer- 
ing Mechanics. Fellow AIAA. 
tScnior Project Engineer. Member AIAA. 


effects of bleed hole size. 1 ’ 14 and location relative to the 
shock.'' 15 Strike and Rippy 1 ' measured the surface pressure 
in the interaction zone of an oblique shock wave impinging 
a turbulent boundary layer over a flat plate, with suction. 
They determined that less suction is required to control sep- 
aration. when applied upstream of the shock. Seebaugh and 
Childs" 1 investigated experimentally the axisymmetric flow in 
the interaction region of the boundary layer inside a duct. 
Contrary to the conclusions of Strike and Rippy.* suction 
within the interaction region was found to be more effective 
in suppressing the effects of separation. 

The purpose of this study is to clarify the effect of bleed loca- 
tion on controlling flow separation in shock wave/boundary- 
layer interactions. Numerical computations are performed to 
determine the flow characteristics in the interaction zone of 
an oblique shock wave/boundary layer, and within a normal 
bleed slot, upstream, downstream, and across the oblique 
shock impingement point. 

Computational Details 

The bleed study is conducted initially for laminar flow, in 
order to allow the fiowfield characterization to be accom- 
plished without the uncertainties associated with turbulence 
modeling. The partial differential equations used to describe 
the flow are the time-dependent compressible Navier-Stokes 
equations in strong conservation law form. The flow is treated 
as a perfect gas with molecular dynamic viscosity given by 
Sutherland's law, and a gas conductivity based on a constant 
Prandtl number. An implicit approximate factorization 
technique 1 ' is used in the numerical solution of the governing 
equations. Prior to the bleed investigation, the code"'- 17 has 
been validated for laminar shock wave/boundary-layer inter- 
actions both for separated and unseparated flows'* by com- 
paring the computational results, with Hakkinen’s experimen- 
tal data.'* 

The solution is obtained in a domain that includes the in- 
teraction region arid the interior of the bleed port as shown 
schematically in Fig. 1. The specified static pressure at the 
bottom of the bleed slot simulates the plenum pressure that 
controls the bleed level for a given slot and flow configuration. 
Zero-order interpolations are applied for the rest of the flow 
variables at the bottom of the bleed slot. The no-slip condi- 
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tion, adiabatic wall, and zero normal pressure gradients are 
specified at the solid surfaces. The impinging shock wave is 
explicitly introduced at the upper boundary using the Ran- 
kine-Hugohiot relations. Zero-order extrapolation of the flow 
variables is enforced at the downstream boundary, which is 
carefully positioned such that all generated or reflected shocks 
exit through it. 

The orthogonal grid spacing was varied in both the stream- 
wise and normal directions to capture the important flow 
characteristics in the shock wave. boundary-layer interaction 
zone and inside the bleed slot. 1 * The computational grid con- 
sisted of (243 x 107) mesh points in the streamwise and 
normal directions, respectively, with grid clustering at the 
plate and around the slot walls. The following grid distribution 
was implemented in the solution of the shock-wave/laminar 
boundary-layer interactions with bleed: 


\Y, = 0.0375 - 

77 </s 107 

(0 

\Y, = 1.125AK 

36.2s./ a 76 

(2) 

AK ( = 1.5625 x 10: 

J 28 ■£ / s 35 

(3) 

1Y, = 1.125AK, , 

1 £ / s 27 

(4) 


where; « 1 represents the bottom, and / = 32 the top of the 
bleed slot. 

The grid inside the bleed slot consisted of (34 x 32) mesh 
points in the streamwise and normal direction, respectively. 
Stretching was used in the streamwise direction, in order to 
cluster the grid points inside, and near the slot as follows: 


IX, = 0.03 / * / ^ 64, / + 63 < / 

IX, = 1.5625 x 10 4 / - 20 ± j s / - 13 

/ + 13 s i s / + 20 

IX, = 1.125*4 M / - 21 at * / - 63 
/ 4 12 * l a / + 1 

&X t = U25A*, / - 12 * is/ 

/ + 21 55 / s / + 62 


(5) 

( 6 ) 

(?) 

( 8 ) 


The flow solution convergence monitoring was based on 
the shear stress at the plate surface r M , as follows: 


r - t!! 


where i - / represents the grid at the center of the bleed 
slot. 
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The specified convergence criterion was e ^ 10 4 for the 
evaluation points at .V < X sh , and e ^ 10 'for the evaluation 
points located at X > to account for the different rates 
for convergence at the different locations from the plate lead- 
ing edge (where X stl is the tested shock impingement point 
location). 

Results and Discussion 

Results are presented for the computed flowfield in the 
interaction zone of an oblique shock/laminar boundary-layer 
interaction with bleed applied through a normal slot at three 
different locations relative to the shock impingement point. 
The freestream conditions are 2.0 Mach number and 2.96 x 
10 5 Reynolds number, based on the shock impingement lo- 
cation, x = 1.96 in., from the flat plate leading edge. The 
shock angle of 32.585 deg corresponds to the test conditions 
of Hakkinen et al.* ,# for the separated flow case. The com- 
putations are performed at the same freestream conditions 
and shock strength, with bleed applied through a normal slot 
that is 0.01 1 in. wide and 0.033 in. deep. The slot width was 
equal to the boundary-layer displacement thickness before 
the interaction, and its depth-to-width ratio of 3 is selected 
based on the experimental results of the optimization study 
by Syberg and Koncsek. 2 The slot was positioned at three 
different locations, upstream, downstream, and across the 
shock impingement point. For across the shock bleed, the slot 
center was positioned directly at the shock impingement point 
with the flat plate surface (4h = 1.96 in.). The downstream 
bleed case was chosen such that the centerline of the bleed 
slot is at the location of the maximum computed surface shear 
stress in the separated flow case without bleed (jr * 2.167 
in ). The upstream bleed location was similarly selected at 
the opposite side relative to the shock impingement point with 
the centerline at x = 1.753 in. from the leading edge. The 
slot base pressure was maintained at the same value (0.42 of 
preshock pressure), which resulted in bleed mass rates of 3.42, 
4.65, and 5.05% of the upstream boundary-layer mass flow 
for upstream, across, and downstream bleed cases, respec- 
tively. Some of the computed results for downstream bleed 
configuration were reported earlier. ■* Here, the results of the 
flowfield computations in the shock/boundary-layer/bleed in- 
teraction zone are presented for the upstream and across the 
shock bleed. The computed results are compared for the three 
bleed locations and for the case with no bleed. 
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Figure 2 shows the pressure distribution near the interaction 
zone for the no bleed case. The incident shock and the sep- 
aration and reattachment shocks are soen in this figure. Figure 
3 shows the plate surface pressure distribution of the three 
bleed cases and for no bleed. The wall pressure distribution 
shows a significant reduction in the length of the interaction 
region compared to the nonbteed case, and an elimination of 
the pressure plateau in all three bleed configurations. Fur- 
thermore, the pressure ratio downstream is slightly increased 
over the nonbleed interaction case. These observations agree 
with previous reports on the effect of bleed, as discussed by 
Hamed and Shang.' Figure 3 clearly shows that there is a big 
difference in the static pressure, at the comer of the flat plate 
surface and the downstream slot wall for the three bleed cases 
This pressure is lowest (p/p. = 1.22) for the upstream bleed 
case, and highest (p/p. = 1.5) tor the downstream bleed case 
with the across the shock pressure value in between (p/p. = 


1.36). As for the static pressure distribution downstream of 
tne bleed slot, there is an abrupt pressure drop for a short 
distance followed by a gradual pressure rise in the upstream 
bleed case. However, in the downstream bleed case, the static 
pressure over the plate surface drops downstream of the bleed 
slot, whereas in the across the shock bleed case the pressure 
rises first, then drops slightly. 

One can see that in both the upstream and across the shock 
bleed cases, the pressure starts to drop before the slot opening 
and reaches a value of p/p. - 0.8 at the intersection of the 
flat plate with the upstream slot wall. In the case of down- 
stream bleed, the pressuie initially builds up on the plate 
surface upstream of the slot opening, reaching a value slightly 
. j « an of the plateau pressure in the case of sepa- 
rated flow without bleed (p/p. - 1.2), followed by a very 
sharp drop to 0.8 at the upstream slot wall comer. In all three 
cases, there is an initial moderate pressure drop over the first 
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part of the slot opening followed by a sharp rise across a shuck 
emanating inside the slot as shown for the upstream bleed 
pressure contours of Fig, 4. 

The Mach number contours are presented in Figs. 5-7. The 
contours indicate a choked flow in all three bleed eases. The 
subsonic boundary-layer flow entering the bleed slot initially 
accelerates due to the separation bubble at the upstream wall 
ot the bleed slot that acts as a converging-type Constriction.. 
After passing through the convergent section of the constric- 
tion. the now choked flow goes through a brief expansion 
region where the flow is further accelerated. The flow then 
decelerates towards the slot exit as it reaches an effectively 
constant area downstream of the separation bubble. 

The velocity distributions across the slot opening are shown 
in Figs. 8- 10. Generally the normal velocity component in- 
creases gradually across the slot opening and reaches a plateau 
. ~~ somewhere j n the downstream half of the bleed 

slot, with a slight increase further near the downstream slot 
wall in the two cases of across the shock and downstream 



Fig. 5 Mach number contours for upstream bleed. 



Hg. 6 Mach number contours for across the shock bleed. 



Fig. 7 Mach number contours for downstream bleed. 
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(Y-Y*)/Q*«0 



x/x* 

Fig. 10 Velocity distribution across slot opening for downstream bleed. 
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X/X„ 

Fig. 1 1 Comparison of the skin friction coefficients for the different 
bleed configurations^ 


bleed. On the other hand, the streamwise velocity component 
reaches a local ulU, of 0.5-0.6 upstream of the slot shock, 
then drops linearly reaching slightly negative values (-0 02 
to -0.1) near the downstream slot corner. Typically, in the 
computations of flowfields involving bleed, the bleed mass 
flux is specified while the tangential velocity is set equal to 
zero across the bleed zone. The results of Figs. 8- 10 suggest 
that the neglect of the streamwise velocity component in the 
bleed zone of shock wave/boundary-layer interaction can con- 
stitute a serious source of error in the computations. This 
might be a major cause of discrepancies between the com- 
puted and experimental results downstream of the interaction 
bleed zone in supersonic inlets. 5 ‘ 7 
Figure 1 1 presents the computed skin friction distribution 
over the flat plate surface for the three bleed locations and 
for the no bleed separated flow case. The reduction in the 
length of the interaction zone d ue to bleed is clear in this 



b) u/u 

Fig. 12 Comparison of the velocity profiles downstream of the in- 
teraction zone for the different bleed configurations. 



X/X„ 


Fig. 13 Comparison of separation region for the different bleed con- 
figurations. 
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4 1 . 



Fig. 14 Comparison of the boundary-layer displacement thickness 
for the different bleed. configurations. 



Fig. 15 Comparison of the boundary-layer momentum thickness for 
the different, bleed configurations. 


figure for all three cases. However, a small flow separation 
zone as indicated by the negative skin friction coefficient can 
be observed for the upstream and downstream bleed cases, 
but not for across the shock bleed. The friction coefficients 
are higher for all three bleed cases downstream of the inter* 
action zone in comparison to the no bleed friction, indicating 
fuller velocity profiles. The highest friction coefficient value 
downstream of the interaction is for the downstream bleed, 
and the lowest for the upstream bleed. These differences in 
the velocity profiles appear to persist for a long distance down- 
stream of the interaction as can be seen in the velocity profiles 
of Fig. 12. 7 y 

Figure 13 compares the extent of the separation regions as 
determined from u - 0. One can observe that the separation 
bubble height is reduced to 32 and 25% of its value without 
bleed for the cases of upstream and downstream bleeds, re- 
spectively. Across the shock bleed completely eliminates sep- 
aration as indicated by the skin friction results of Figs. 11 and 
12. The boundary layer's displacement and momentum thick- 
ness are presented for the three bleed locations and for no 
bleed in Figs. 14 and 15. One can see that the largest reduction 
in 6* and 6 are realized in the downstream bleed case and the 
least reduction for upstream bleed. 


Conclusions 

A numerical investigation was conducted to determine the 
effect of bleed location on the flow in an oblique shock/bound- 
ary-laycr interaction. The implicit factored scheme of Beam 
and Warming 15 was used to obtain numerical solutions of the 


two-dimensional compressible Navier-Stokes equations in 
strong conservation law form . The investigation was restricted 
to laminar flows to avoid the uncertainties associated With 
turbulence modeling. The solution domain extended inside 
the normal bleed slot, and computational grid was clustered 
near the surfaces. Three bleed locations upstream, across, 
and downstream of the shock impingement point were in- 
vestigated for the same slot size and static pressure at the 
bottom of the bleed port. Results were presented for pressure 
and Mach number throughout the interaction region and in- 
side the bleed port, and for the surface pressure and skin 
friction distribution. According to the computed results, the 
across-shock bleed proved to be the most effective in pre- 
venting separation and shortening the interaction length, 
whereas the upstream bleed slot proved to be the least ef- 
fective. The results indicate that the downstream bleed slot 
location produced the most favorable velocity profiles and 
skin friction recovery downstream of the interaction region, 
as well as the maximum reduction in boundary-layer mo- 
mentum and displacement thickness. The results suggest that 
a combination of across-shock and downstream bleed slots 
might provide the two advantages, namely no flow separation 
in the interaction region and fuller boundary-layer profiles 
downstream of the interaction zone. 
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Investigation of Shock/Turbulent Boundary-Layer 

Bleed Interactions 


A. Hamed,* S. Shih,t and J. J. Yeuant 
University of Cincinnati, Cincinnati, Ohio 45221 


A numerical investigation was conducted to determine the effect of bleed on oblique shock wave/turbulent 
boundary-layer interactions (SWBLI). The numerical solutions to the compressible Navier-Stokes equations 
reveal the flow details throughout the interaction zone and inside the normal bleed slot. Results are presented 
for an incident oblique shock of sufficient strength to cause boundarv-layer separation on a flat plate in the 
absence of bleed at a frees tr earn Mach number of 2.96 and Reynolds number of 1.2 x Witt, with bleed applied 
across the shock impingement location over a range of bleed mass flow rates corresponding to different values 
of plenum pressures. The results indicate a complex flow structure with large variations in both the normal 
and tangential flow velocities across the bleed slot. The flow entrainment into the slot is accompanied by an 
expansion-compression wave system with a bow shock originating inside the bleed slot. Increasing the bleed 
mass flow by decreasing the plenum pressure caused an initial decrease, then a later increase in boundary- 
layer momentum and displacement thickness downstream of the interaction. 


Introduction 

T HE control of flow separation in shock boundary-layer 
interactions can be accomplished through flow suction 
(bleed) in the interaction zone. In the case of mixed compres- 
sion supersonic inlets, the bleed system design is Critical to 
the efficient and stable operation of the system. Hamed and 
Shang 1 reviewed the existing experimental data for shock wave 
boundary-layer interactions in supersonic inlets and other re- 
lated configurations. According to this survey, there is enough 
experimental evidence 2-7 to indicate that local bleed can con-, 
trol flow separation in shock wave/boundary-layer interac- 
tions. There are disagreements, 1 however, among the differ- 
ent experimental results regarding the effects of bleed location 
relative to the shock 2-4 - 7 and of bleed hole size. 7 -* 

Strike and Rippy 2 obtained surface pressure measurements 
with porous bleed applied across and upstream of the shock 
impingement. Based on the minimum mass flow required for 
the static pressure downstream of the interaction to reach the 
theoretical pressure ratio, they determined that less suction 
is required when applied upstream. Seebough and Childs 3 who 
conducted their oblique shock wave boundary-layer bleed study 
over the inner surface of a cylinder, with a cone as a shock 
generator, concluded that across the shock bleed was more 
effective in controlling boundary-layer separation. Hingst and 
Tanji 4 and Benhachmi et al. J used pressure taps to measure 
the surface pressure, and hot wire to measure the bleed mass 
flow distribution through normal holes in SWBLI. Their re- 
sults indicated that the bleed mass flow distribution followed 
the surface pressure distribution. In addition, Hingst andTangj 4 
used pitot probe surveys and surface flow visualization to 
compare the performance of bleed locations upstream, across, 
and downstream the shock. According to this study, across 
the shock and downstream bleed produced similar results in 
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terms of removing flow separation, shortening the interaction 
length, increasing the pressure gradient, and producing fuller 
velocity profiles downstream of the interaction. Lee and 
Leblanc* presented their experimental results from taps, schli- 
eren, pitot, and static probes for two levels of suction applied 
through porous wall sections across the shock impingement 
point. The poor performance of the weak suction was attrib- 
uted to the effects of surface roughness. 

Fukuda et al. 7 determined experimentally the change in 
boundary-layer characteristics across oblique shock wave in- 
teractions with bleed over the centerbody and cowl of a super- 
sonic inlet. This study was conducted at different bleed mass 
flow rates, bleed hole 3izes, bleed locations relative to the 
shock, and bleed hole inclination relative to the surface. Based 
on the change in boundary-layer displacement and momen- 
tum thicknesses, and in the transformed form factor for the 
different bleed configurations, they concluded that bleeding 
upstream and downstream of the interaction was preferable 
to bleeding across the shock. These results were unchanged 
by the normal bleed hole size which was found to have a 
negligible effect on the results. 

Most viscous flow computations in supersonic inlets 9- 13 re- 
quire some knowledge of bleed flow rates and/or pressure 
distribution in the bleed regions. Comparisons of internal flow, 
computational results with the experimental measurements in 
supersonic inlets revealed reasonable agreement between the 
computed and measured surface pressures upstream of the 
ramp bleed. However, discrepancies in the predicted shock 
locations and velocity profiles downstream of shock bound- 
ary-layer interactions with bleed suggest deficiencies in the 
bleed modeling. Other investigators actually solved the flow 
through the boundary layer and into the bleed slot, providing 
insight into the bleed flow region;'*- 17 This approach elimi- 
nates the need for experimentally derived bleed models for 
supersonic inlet flow simulations, but can be very costly com- 
putationally. 

Recently, some research efforts have been directed at im- 
proving bleed models. Paynter et al. 1 * used the rough wall 
algebraic turbulence model of Cebcci and Chang 1 * to simulate 
the overall effect of bleed on the growth of the boundary 
layer. The proper values of the roughness parameters were 
determined through matching the computational and test re- 
sults for near wall velocity distributions downstream of the 
bleed band for seven bleed configurations, with shock bound- 
ary-layer interactions in some cases. The results of this in- 
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vestigation demonstrate a strong dependence of the surface 
roughness On the fraction of the upstream boundary-layer 
mass flux removed. High roughness values were found for 
bleed rates between 3-15% of the upstream boundary-layer 
mass flux, but were insignificant for lower and higher bleed 
rates. Choked bleed flow through slanted holes gave the best 
boundary-layer profile in terms of flow separation control. 

Hamed and Lehnig 20 - 21 conducted a numerical study in an 
incident oblique shock/laminar boundary-layer interaction, with 
bleed through a normal slot. By obtaining the flow solution 
in a domain that included the shock boundary-layer and ex- 
tended inside the bleed slot, they were able to reveal the 
various regimes in this complex flOwfield. Their results dem- 
onstrated a large variation in both the normal and tangential 
flow velocity distribution across the bleed slot opening at the 
plate surface. A separated flow region inside the bleed slot 
in the computed results suggests a reduction in the flow coef- 
ficient and total pressure. Hamed and Lehnig 21 presented the 
flowfields for three locations of the normal bleed slot, up- 
stream, across, and downstream of the shock impingement 
point. The results demonstrated that for a given plenum pres- 
sure at the slot exit, the shock-induced flow separation on the 
plate was eliminated in the across the shock bleed configu- 
ration and reduced in the upstream and downstream bleed 
cases. Edwards and McRae 22 also computed the flowfield for 
an oblique incident shock, on a laminar boundary layer, with 
bleed through a normal slot across the shock impingement 
point. They compared the computed surface pressure distri- 
butions for a specified uniform normal bleed velocity across 
the slot to the plenum case, and demonstrated that the first 
case underpredicts the initial flow expansion into the slot. The 
solution domain for the plenum case was similar to that of 
Hamed and Lehnig, 20 - 21 but a uniform suction velocity was 
specified at the slot exit. 

The discrepancies among the different experimental studies 
is an indication of the complexity of the flow in these config- 
urations. Based on the comparison of their computational 
results with the experimental data of supersonic inlet flow- 
fields, Reddy et al.' 3 stressed the need for a detailed study of 
the effect of the individual bleed ports. Bleed optimization 
can only be accomplished through a parametric investigation 
in which the bleed conditions are changed systematically. The 
large number of parameters and the difficulties in obtaining 
accurate flow measurements in the interaction zone precludes 
a complete experimental investigation. 

The purpose of the present study is to conduct a numerical 
investigation to characterize the flowfield in oblique shock 
wave/turbulent boundary-layer interactions with bleed, and 
to provide a basic understanding of the bleed flow controlling 
mechanisms. In the investigated configuration, flow suction 


(bleed) is applied through a normal slot located across the 
incident oblique shock on a flat plate turbulent boundary 
layer. The region inside the slot from the bleed surface to the 
plenum exit was included in the solution domain in order to 
realistically model the bleed flow and its variation with the 
plenum pressure. The results of the flow computations with 
different bleed mass flow rates reveal in detail the various 
flow regimes in the interaction region and within the bleed 
slot. 

Computational Details 

The PARC code 23 was used in the present investigation 
after conducting validation studies for shock wave/turbulent 
boundary-layer interactions, 24 with and without flow sepa- 
ration. The performance of four different codes with different 
turbulence models was assessed in terms of the agreement of 
their computed results with the experimental data of Refs. 25 
and 26. The required grid refinement and computer time to 
reach this agreement were compared for incident shock on a 
flat plate boundary layer, and for flow over a compression 
comer. The prediction of the separation length, the surface 
pressure propagation upstream of the interaction, and the 
recovery of the skin friction downstream of the interaction 
received special attention. From the four assessed codes, the 
two 23 27 that use the implicit approximate factorization tech- 
nique of Beam and Warming 28 to advance the solution, per- 
formed best both in terms of computer time and grid inde- 
pendence. Visbal’s code 27 was the most robust, required the 
least number of time steps to reach convergence, and its com- 
puted results were consistently better than those obtained 
from the other codes with algebraic turbulence model. How- 
ever, Awa et al. 29 proved the two equation k-e model to be 
superior to the algebraic models for applications with sepa- 
rated flows. Consequently, the PARC code 23 was used in the 
present study because among the codes with the two equation 
turbulence models, it performed best, in terms of agreement 
with the experimental results and required the least grid re- 
finement to achieve this agreement. 24 In addition, the PARC 
code was validated 30 through comparisons with experimental 
data for several flow configurations including supersonic flow 
over compression with separation, and has been used in sev- 
eral supersonic and hypersonic inlet flow simulations. 12 - 14 - 31 

Central differencing explicit and implicit dissipation schemes 32 
are used in the PARC code for the solution of the Reynolds- 
averaged Navier-Stokes equations in strong conservation form. 
The Beam and Warming factorization algorithm coupled with 
Pulliam’s pentadiagonal formulation form the basis of its im- 
plicit LU-ADI style solution scheme. The turbulence model 
used in the PARC code is based on Chien’s low Reynolds 
number k-e model 33 with Nichols 30 modifications to add com- 
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pressibility effects. For the simulation, flow was considered 
to be turbulent throughout the calculation domain ; no attempt 
was made to model transition and/or relaminarizution. 

Refenring to Fig. 1, the solution domain used in the two- 
dimensional flow simulations extends upstream of the flat 
plate leading edge, downstream of the Shock wave/turbulent 
boundary-layer interaction region, and inside the bleed slot, 
where the specified plenum pressure at the bottom controlled 
the amount of bleed mass flow. The incident oblique shock 
crosses the upper boundary, all other wave systems including 
the reflected, and any separation and reattachment shocks, 
cross the outflow boundary. Uniform freestreum flow con- 
ditions are applied at the inflow boundary and all variables 
were extrapolated at the outflow boundary. No slip adiabatic 
flow Conditions are applied at the plate and bleed walls. The 
static pressure is specified and first Order extrapolation is 
applied for the rest of the flow variables at the bottom of the 
bleed slot. The location of the incident shock is fixed at the 
upper boundary with freestream and postshock conditions 
specified upstream and. downstream. The slot was centered 


. around the shock impingement location 1 ft from the flat plate 
leading edge. The grid spacing was varied in both x and y 
directions with clustering around both bleed walls and at the 
plate surface as shown in Fig. 2. The computational results 
were obtained with Ay min = Ax mm = 0.25794 x 10 - 4 ft. 
corresponding to y * * 2.0 at x - 0.9 ft, using a 308/68 grid 
over the plate surface and a 75/68 grid inside the bleed slot, 
based on the grid refinement study results of Ref. 24. The 
computations for a typical bleed case required 5000 local time 
steps at 0.3 Courant-FriedriChs-Lewy (GFL) number to reach 
steady-state solution based on six orders of magnitude re- 
ductions in the averaged rms error in the flux. 

Results and Discussion 

The Computations were performed at the incoming flow 
conditions of M, = 2.96, Re, = 1.2 x lO’/ft, T„ = 437°R 
and P„ = 54.7 psi, and an impinging oblique shock angle of 
25.84 deg. These are the same as the experimental test con- 
ditions of Law 23 for the separated flow case at a deflection 
angle 8 of 7.93 deg. In his experimental study. Law used oil 
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Fig. 2 Computational grid in the interaction region, bleed width a 0.0213 ft. 
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Fig. 3 Surface pretture distribution in the interaction region without 
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Fig. 4 Friction coefficient distribution in the interaction region with* 
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flow visualization to reveal non-two-dimensional effects of the 
corner on the flat plate interior region. He measured the 
separation length for different shock generation spans and 
used these results to determine the proper reduced span to 
minimize the influence of the side wall interactions on the flat 
plate centerline. Figures 3 and 4 compare the computed sur- 
face pressure and skin friction without bleed to Law’s cen- 
erline experimental data. One can see that the predicted 
separation length and surface pressure distribution are in very 
good agreement with the experimental results. The compu- 
tations accurately predict the extent of the upstream influence 
of the shock boundary-layer interactions on the surface pres- 
sure and slightly overpredict the separation length. 

The computations were performed under the same flow 
conditions with a bleed slot normal to the plate at different 
plenum pressure ratios (0.95 s p b lp x s 1.95). For an inviscid 
pressure ratio across the reflected shock p 2 lp* equals 2.962, 
and based on the downstream pressure p 2 the bleed plenum 
pressures are set to vary from choked to unchoked slot con- 
ditions (0.321 < p b lp 2 < 0.658). The bleed slot width was 
0.0213 ft, which is equal to 1.617 times the boundary-layer 
thickness upstream of the interaction at x = 0.9 ft. The details 
of the ftowfleld computational results near the bleed slot en- 
trance are presented in Figs. 5-7 for a plenum pressure ratio 
of 1.5. They indicate a complex flow structure in which an 
expansion-compression system is associated with the flow en- 
trainment into the slot opening. The initial flow turning is 
weak, then the flow is turned more strongly, in the second 
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half before crossing the bow shock that originates inside the 
slot. The Mach number contours show a stagnation point 
behind the bow shock near the plate comer, followed by an 
expansion fan. The flow inside the bleed slot is mostly con* 
fined to a strip adjacent to the back wall of the bleed slot, 
where the flow velocities reach supersonic values. A flow 
recirculation zone much larger than the laminar case 21 oc- 
cupies most of the bleed slot. The computational results at 
different pjenum pressure settings show that the dual sepa- 


ration and reattachment shock system is replaced by a single 
reflected shock structure, indicating no flow separation, for 
plenum pressure ratios sl.S. The corresponding variations 
in static and total pressures, static temperature, and Mach 
number across the bleed slot are presented in Figs. 8-11. 
These figures show the sudden change in the flow properties 
across the bow shock. One can observe a reduction in both 
the static and total pressure over the plate surface upstream 
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of the slot opening due to the flow entrainment into the slot. 
Downstream of the slot opening, both the static and total 
pressure increase more sharply compared to the no bleed case. 

Figure 12 presents the velocity distribution across of the 
bleed slot opening for different bleed flow rates. Large var- 
iations in both the normal and tangential velocity components 
can be seen across the slot opening with some mass injection 
near the upstream corner at the high-pressure plenum set- 
tings. The tangential velocity reaches values as high as 90% 
of the freestream velocity before dropping sharply behind the 
bow shock near the back wall. 

The variations in the static pressure at the bleed slot exit 
with the calculated bleed mass flow are presented in Fig. 13. 
The corresponding changes in the mass averaged total pres- 
sure at the bleed slot inlet and exit are presented in Fig. 14. 
In these figures, the bleed mass flow is normalized with re- 
spect to the mass flow in the incoming boundary layer at 
x = 0.9, and the static and total pressures are normalized 
with respect to the freestream values. With the static pressure 
nearly constant across the bleed slot opening, except for the 



Fig. 16 Near wall velocity profiles, downstream of the interaction at different bleed pressure ratios. 
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Fig. 17 Velocity profiles downstream of the Inte raction at different Meed p resto re ratios. 
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sudden rise across the bow shock, the total pressure increase 
across the bleed slot opening is mainly associated with the 
increased flow velocities as the bleed mass flow increases. 
However, at the slot exit, the mass averaged total pressure 
is mostly influenced by the value of the static pressure due to 
the greater effects of the flow recirculation before choking. 
The corresponding variation in the boundary-layer displace- 
ment and momentum thickness are presented in Fig. 15. The 
results indicate that the boundary-layer displacement and mo- 
mentum thicknesses downstream 6f the interaction initially 
decrease, then slightly increase with the increase in the bleed 
mass flow. According to Fig. 15, the maximum reduction in 
the boundary-layer displacement and momentum thickness 
downstream of the interaction are achieved at plenum pres- 
sure ratios of 1.2 and 1.5, respectively. The corresponding 
bleed mass flows are 27.38 and 24.51% of the inco ming 
boundary-layer mass flow upstream of the interaction. 

■Hie near wall velocity profiles downstream of the inter- 
action are presented in Fig. 16 for three cases with plenum 
pressure ratios of 1.2, 1.5, and 1.75. This figure shows that 
the velocities near the wall are higher for the case p/p. = 1.5 
than for both higher and lower bleed mass flow cases, p/p. 
= 1.2 and p/p. = 1.75. This can be attributed to the balance 
between the increases in the flow velocities and bow shock 
strength with the increased bleed mass flow. Figure 17 shows 
the variation in the flow velocity profiles downstream of the 
interaction Outside the boundary layer for the three bleed 
cases corresponding to plenum pressure ratios of 1.2, 1.5, and 
1.75 and for the no bleed case. This figure illustrates the 
reflected shock structure, with the strong bow shock closest 
to the wall, followed by the expansion fan, then the weaker 
separation shock further out. Figure 17 shows that as the bleed 
mass flow increases, the reflected shock system moves closer 
to the wall. 


Conclusions 

The objective of this work was to provide a basic under- 
standing of the flow in the shock boundary-layer interactions 
and of the mechanisms for its control. Numerical computa- 
tions were conducted in oblique shock wave/turbulent bound- 
ary-layer interaction to reveal in details the flow character- 
istics, and how they are altered by bleed. Flow suction (bleed) 
was accomplished through a normal slot across the shock 
incident point. The bleed results reveal a complex flow struc- 
ture in the interaction zone and inside the bleed slot. The 
flow accelerates across the top of the bleed port, crosses a 
bow shock formed near the downstream comer, then goes 
through an expansion fan at the comer. A large recirculation 
region dominates the flow inside the bleed slot, and the bleed 
mass flow is limited to a narrow region adjacent to the down- 
stream bleed wall, where the velocity reaches supersonic val- 
ues. The effect Of the bleed mass flow rate was investigated 
by changing the specified plenum pressure at the bottom of 
the slot and comparing the boundary-layer characteristics 
downstream of the interaction. As the bleed mass flow in- 
creased, the boundary-layer momentum and displacement 
thickness initially decreased, then slightly increased down- 
stream of the interaction. In terms of the boundary-layer char- 
acteristics downstream of the interaction, the best perfor- 
mance was obtained when 24-27% of the boundary-layer 
ma«« °. 'vi spctream of the interaction was removed. 
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Abstract 

A numerical investigation was conducted to study the 
effect of bleed configuration on oblique-shock wave/turbulent 
boundary-layer interactions. Bleed is applied through a 
normal slot across the shock impingement location. The 
numerical solution to the compressible Navier-Stokes 
equations is obtained for the turbulent flow throughout the 
interaction zone and inside the bleed slot for bleed mass flow 
rates up to 57% of the boundary layer mass flow upstream of . 
the interaction. The results indicate that the bleed slot 
performance improves as the slot width decreases and the 
length to width ratio increases. This is reflected as an . 
increase in the bleed discharge coefficient and total pressure, 
and a reduction in the boundary layer momentum and 
displacement thickness downstream. 

Introduction 

The control of shock wave boundary layer interaction . 
through bleed (suction) is essential for the efficient and stable 
operation of supersonic inlets. Bleed . educes the unfavorable 
effects associated with these interactions such as boundary 
layer thickening and flow distortion, and is especially critical 
for interactions with shock strong enough to cause separation. 
Hamed and Shang [1] reviewed the existing data base on 
shock wave boundary layer interactions in supersonic inlets. 

In spite of the very limited number of investigations with . 
local flow measurements in the interaction region, their 
conclusions regarding the effect of bleed configuration were 
not in agreement. Fukuda et al. [2] reported no significant 
change in the bleed system performance with the change in 
bleed hole size. Wong [3] on the other hand reported 
significant hole size effects and recommended hole sizes 
below the boundary layer thickness for better bleed 
performance. Similarly, contradictions can be found 
regarding the effect of bleed location relative to the shock [4, 

5]. 

Different approaches have been tried to model bleed 
without resolving the flow through the individual bleed holes. 
Ccbcci's algebraic turbulence models [6, 7] for rough wall 

boundary layers, have been used by Abrahmson et al. [8], 

* 
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Benhachmi et al. [9] and Paynter et al. [10] to simulate the 
bleed hole roughness effects in the numerical solutions to the 
Navier-Stokes equations. Empirical boundary conditions for 
the mass flux distribution across the bleed zone in the shock 
boundary layer interactions were used with this turbulence 
model to simulate the bleed influence in the numerical 
computations of references [8] and [9], Paynter et al. [10] 
used the experimental data for different bleed configurations 
to determine the appropriate value of roughness [7] that 
models the overall effect of bleed on the boundary layer 
growth rate. They determined that the roughness parameter 
depends on the fraction of the upstream boundary layer mass 
flow removed. 

Hamed et al. [11-13] and Rimlinger ct al. [14] 
conducted numerical studies of shock wavc/boundary 
layer/bleed interactions which simulated the flow field 
throughout the interaction zone and inside the bleed port. 
These studies revealed a complex pattern of expansion and 
compression waves across the bleed opening for selected slot 
[11-13] and hole [14] configurations. Both Hamed ct al. [12] 
and Rimlinger et al. [14] studied the effect of bleed location 
relative to the shock on the resulting flow field. In addition, 
Hamed et al. [13] predicted the variations in the boundary 
layer characteristics downstream with the bleed mass flow. . 

A complete bleed data base is critical to achieving inlet 
bleed system design objectives of eliminating or reducing the 
effects of flow separation in SWBLI, while minimizing bleed 
drag and bleed mass flow. A large number of parameters 
such as the local Mach number and shock strength, the bleed 
hole size, slant angle and the bleed location relative to the 
shock can influence the performance. The various bleed 
studies addressed different aspects of this Complex 
phenomena and presented results for one or more of the bleed 
performance parameters. The test results of Syberg and 
Koncsek [15] consisted of the measured bleed mass flow 
variation with plenum pressures and the bleed drag for 
choked holes. The results were presented for different bleed 
hole diameters, aspect ratios and slanting angles at free, 
stream Mach numbers ranging between 0.8 and 2.19. The 
test cases did not include shock interactions or pressure 
gradient. Fukuda et al. [2] obtained boundary layer 
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measurements downstream of boundary layer interactions on 
the center body and cowl of a supersonic inlet. Over the 
range of the tested bleed mass flows, the transformed form 
factor downstream of the interaction initially decreased then 
reached a plateau or slightly increased as the bleed mass flow 
increased. As expected, the slanted holes produced the lowest 
form factor but the hole size effect was insignificant. Hingst 
and Tangi [16] studied bleed configurations consisting of 
normal bleed hole rows in regions upstream, across and 
downstream of the impinging shock. Their results indicated 
that fuller boundary layer profiles, shorter interaction length . 
and higher pressure gradients were obtained with the bleed . 
region across the shock. They also used hot wire to measure - - 
the bleed mass flux distribution through the porous wall 
which was found to change in a pattern similar to the surface 
pressure. - 

Based on the previous discussion, the data base required . 
to optimize bleed system design is not complete and does not 
provide sufficiently detailed information. The purpose of the 
present numerical study is to investigate the effect of bleed 
configuration on the flow throughout the shock boundary 
layer interaction and inside the bleed port. A parametric 
study was conducted in which bleed configuration parameters 
such as the bleed slot width and bleed aspect ratio Were 
changed systematically to determine their effects over a wide 
range of bleed plenum pressures. The results are presented . _ 
for the flow field and also for the bleed performance 
parameters in terms of the discharge coefficient, total 
pressure recovery and boundary layer characteristics 
downstream. 

Comnutational Details 

The PARC code [17] with a two-equation k-e turbulence - 
model was used in this bleed study. The turbulence model is 
based on Chien's low Reynolds number k-e model [18] with 
Nichols [19] modifications to add compressibility effects. 

This code has been used in several supersonic and hypersonic 
inlet flow computations [20-2 1]. 

Numerical flow simulations were conducted for bleed 
applied through a slot normal to the wall in the interaction 
region of an externally generated incident oblique shock and 
a turbulent boundary layer. The flow was considered to be 
turbulent throughout the calculation domain; no attempt was 
made to model transition and/or relaminarization. As in our 
previous studies [11-13], the solution domain used in the two 
dimensional flow simulations extends upstream of the flat 
plate leading edge, downstream of the interaction region and 
inside the bleed slot, where the specified plenum pressure at 
the bottom controlled the amount of bleed mass flow. The 
incident oblique shock crosses the upper boundary, and all 
other wave systems including the reflected, and any 
separation and reattachment shocks, cross the outflow 
boundary. Uniform freestream flow conditions are applied at 


the inflow boundary and all variables are extrapolated at the 
outflow boundary. No slip adiabatic flow conditions arc 
applied at the plate and bleed walls. The static pressure is 
specified and first order extrapolation is applied for the rest of 
the flow variables at the bottom of the bleed slot. The 
location of the incident shock is fixed at the upper boundary 
with free stream and post shock conditions specified upstream 
and downstream. 

A 308x68 grid was used over the plate surface with 
Aymin = 0.25794 x 10* 4 ft. corresponding to y + =2.0 
upstream of the interaction at x - 0.9 ft. An 89x68 grid was 
used inside the slot which was centered around the shock 
impingement location at 1 ft from the flat plate leading edge. 
The grid spacing was varied in both x and y directions with 
clustering around both bleed walls and at the plate surface 
[13]. The computations for a typical bleed case required 5000 
-local time steps at 0.3 CFL number to reach steady state 
solution based on six orders of magnitude reductions in the 
averaged root mean square error in the flux. 

Results and Discussion 

The computations were conducted at the incoming flow 
conditions of = 2.96, Re^ = 1.2 x l0 7 /ft, and impinging 

oblique shock angle of 25.84° (wedge angle of 7.93°) 
corresponding to the test conditions of reference [22] for the 
separated flow case with no bleed. The flow characteristics 
“inside the normal bleed slot, and throughout the interaction 
were determined for different values of bleed mass flow, slot 
width, D, and length L. The pressure at the bottom of the 
bleed slot, p^, was varied between 1.95 and 0.42 times the 
free stream values to obtain normalized bleed mass flows 
ranging from zero up to 0.5679 of the incoming boundary 
layer mass flow rate. The results are presented and compared 
for three different bleed slot width values of 0.0106, 0.0213 
and 0.0426 ft. which represent 0.8085, 1.617 and 3.234 times 
the boundary layer thickness upstream of the interaction, and 
for bleed slot length to width ratios of 3.05 and 6.10. 

Typical results of the flow field computations are 
presented in Fig. 1 for two bleed cases with slot widths of D/S 
* 0.8085 and 1.617 at pjj/pi oc =0.845. Inside the slot, a 
significant part of the width is occupied by a large 
recirculaton zone while the bleed flow near the back wall 
reaches supersonic velocities. The pressure contours indicate 
a complex wave system across the bleed slot opening, starting 
with an expansion fan at the upstream slot comer, then a bow 
shock originating inside the slot. A separated flow region 
over the plate surface upstream of the narrow slot is indicated 
by the Mach number contours and the island of high pressure . 
contours above the slot at this pressure level. This effect is 
greatly reduced in the wider slot with the higher bleed mass 
flow, where the flow reaches higher Mach numbers upstream 
of the outgoing slot shock. Figure d presents the Mach 
number contours for the narrow slot at different bleed plenum 
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pressures while Fig. 3 presents the corresponding variation in 
the flow properties across the slot opening. One can see that 
as the plenum pressure is reduced and the flow expands to 
higher Mach numbers at the slot opening, the shock 
originating in the slot becomes stronger. Eventually, the flow 
conditions across the slot opening remain unchanged with 
further plenum pressure reductions below Pb/Ploc^- 383 1 . 
Figure 2 indicates that at these low plenum pressures, the 
flow separates along the slots downstream wall, effectively 
producing a converging flow path leading to supersonic flow 
diffusion inside the slot from M = 2.3 to l.S. On the other 
hand, at the highest plenum pressure, lVPloc =1 - 098 > Fi S- 3 
indicates that no shock originates inside the slot, but that the 
mass flux out of the slot opening behind the upstream comer 
causes a weak oblique shock to form there. . 

The detailed flow field results for the different bleed . 
configurations were further processed to characterize the 
bleed system performance in terms of integrated parameters 
such as the discharge coefficient, total pressure recovery and 
the boundary layer displacement and momentum thicknesses 
downstream of the interaction. The bleed mass flow for the 
individual slot is characterized in terms of the discharge 
coefficient,, which is defined as follows: 
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In the above definition, the bleed mass flow is normalized 
with respect to the ideal mass flow if isentropic flow was to 
expand to sonic conditions from a total temperature equal to 
that of the free stream and a total pressure equal to the local 
static pressure. With the bleed slot centered at the shock 
incident point at the surface, the post shock static pressure, 
Ploc> ' vas used in the bleed discharge coefficient calculations . _ 
At the incident shock wedge angle of 7.93° and an incoming 
free stream Mach number of 2.96, the post shock static and 
total pressure ratios are Pioc^Pinf = 1-7752 and Pt,loc/Pt,inf = 
6.981 respectively, while the corresponding values after the 
reflected shock are P2^Pinf = 2 96 an( l Pt2^Pt,inf = 

Effect of Slot Width 

Figures 4a and 4b present the variation in the bleed 
discharge coefficient and the normalized bleed mass flow 
with the bleed plenum pressure for different slot widths. 
According to this figure, the discharge coefficient is higher 
for the narrower slots, and the bleed mass flow continues to 
increase for Pt/Pioc < 0.528, reaching a maximum in the 
narrow and intermediate slot at pj/pio^.2. Figure 5 
presents the variation of the mass averaged total pressure at 
the bleed slot inlet and exit with the bleed mass flow for the 
different bleed widths. According to Fig. 5, the mass 
averaged total pressure at the slot opening changes only with 


the bleed mass flow and is independent of the bleed geometry. 

It increases with the bleed mass flow for values above 12% of 
the incoming boundary layer. The decrease in the slot inlet 
total pressure With the bleed mass flow below this value is 
associated with the formation of the slot bow shock at p^/pioc 
< 1.09 and its location relative to the slot back corner as was 
seen in Fig. 3. The behavior of the mass averaged total 
pressure at the slot exit is more complex (Fig. 5). Initially it 
decreases with the bleed mass flow for the narrow and 
intermediate slot, while it remains unchanged for the wide 
slot. Then the total pressure at the slot exit increases for all 
slot sizes and peaks at a different bleed mass flow in each 
case corresponding to Pb/piocS 0.528 (see Fig. 4). The 
subsequent drop in the exit total pressure in Fig. 5 was found 
to be associated with flow separation over slot back wall at 
the lower plenum pressures (see Fig. 2). 

The results for the separated flow length over the plate 
surface, and the boundary layer displacement, and momentum 
thickness downstream of the interaction are presented in Figs. 
6a through 6c. According to these results, both the 
displacement and momentum thicknesses initially decrease 
then start to increase with increased bleed mass flow even 
while the flow separation is reduced or eliminated. The 
narrow slot is found to perform best in terms of the maximum 
reduction in 6 and 5* with minimum bleed mass flow, 

mb / nibi <0.10. To explain this phenomena. Figs. 7a and 
7b present the total pressure distribution across the opening 
and the skin friction variation over the plate surface for the 
three slot widths at approximately, the same bleed mass flow, 

mb /nibi =0.12. According to this figure, the large 
reductions in 6, 8* are closely associated with the weaker 
shock inside the narrow slot, and not with the elimination of 
flow separation over the plate surface. 

Effect of Slot Aspect Ratio 

The pressure, and Mach number distribution, across the 
slot opening for aspect ratios of 3.05 and 6.10 are presented 
in Fig. 8. According to these figures, lower static pressures, 
and higher flow velocities are associated with higher aspect 
ratio. The effect of the slot length to width ratio on bleed 
performance is presented in Fig. 9. The increase in L/D 
increases the discharge coefficient for both the narrow and 
intermediate slots (Fig. 9a) On the other hand, the bleed exit 
total pressure decreases with the increased aspect ratio as 
shown in Fig. 9b. The boundary layer momentum thickness 
downstream of the interaction decreases with L/D in the case 
of the intermediate slot and slightly increases for the narrow 
slot, but is generally not very sensitive to changes in L/D. 

Conclusions 

A numerical investigation was conducted to reveal the 
flow characteristics, and how they are affected by the bleed 
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configuration in an oblique shock wave/turbulent boundary 
layer interaction. Flow suction (bleed) was accomplished 
through a normal slot across the shock incident point. The 
bleed results reveal a complex Wave pattern across the slot 
opening consisting of an expansion wave at the upstream slot 
corner, a bow shock wave originating inside the slot. The _ 
bleed flow reaches supersonic velocities inside the slot with . 
one or two large recirculation regions at the walls depending 
on the plenum pressure. The effect of bleed geometry was 
investigated by changing the slot width and length 
systematically and determining the corresponding flow field 
at different plenum pressures. The narrow slot (D/8 = 0.808) 
gives the best bleed performance with the highest discharge 
coefficient and the lowest momentum thickness downstream 
for the least bleed mass flow. Most importantly, eliminating 
separation on the plate surface was not required to optimize 
the boundary layer characteristics downstream, which was 
achieved at narrow slot bleed mass flow of mt/mjjL <0.1. As 
the slot aspect ratio increases, the bleed discharge coefficient 
increases, but the bleed flow total pressure decreases. The 
boundary layer momentum thickness downstream of the 
interaction was not very sensitive to changes in L/D. 
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Fig. lb. Mach Number Contours, D/6=0.8085 
L/D = 3.05 and P^/p^ = 0.676. 
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Fig. Id. Mach Number Contours, D/6=1.617 
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Fig. 4. Effect of Slot Width on Bleed Discharge Coefficient and Maas Flow Rate, L/D - 3.05. 
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Fig. 5. Effect of Slot Width on Bleed Mass Averaged Total Pre ssure. UD - 3.05 
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Fig. 6. Effect of Slot Width on Separation Length, Displacement and Momentum Thickness Downstream of 
Intetaction at % » 1.13 ft, L/D - 3.0S 
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ABSTRACT 

Hie flow field characteristics are simulated 
numerically in an oblique shook wave/turbulent boundary 
layer interactions with six different bleed slot 
configurations. The strong conservation-law form of the 
two-dimensional compressible Navier-Stokea equations and 
the k-€ equations are solved throughout the interaction 
region and inside the bleed slot The computed results are 
presented for a normal and 20* slanted bleed slots at three 
different locations, upstream, across and downstream of the 
impingement point of an oblique shock of sufficient 
strength to cause boundary layer separation without bleed. 
The detailed flow characteristics in the interaction zone and 
inside the bleed slot are compared for the different bleed 
slot configurations. The resulting surface pressure and 
shear stress distributions as well as the boundary layer 
characteristics downstream of the interaction region are 
also presented for the six bleed configurations at different 
bleed mass flows up to choking. 




Air bleed systems are used for controlling the 
shock wave boundary layer interactions when operating at 
supersonic speeds. Proper bleed system design is 
particularly important for the efficient and stable operetion 
of mixed compression supersonic inlets [1]. The 
fu nd a men tal objectives of supersonic inlet bleed system 
design are to provide good aerodynamic flow 
characteristics with minimum boundary layer bleed [2, 3], 

Comparisons of internal flow computational 
results [S, 6, 7] with the experimental measurements in 
supersonic inlets [7, 8] revealed reasonable a g reement 
between the computed end measured surface pressures 
upstream of the lamp bleed. However, discrepancies in the 
predicted shock locations end velocity profiles were 
observed downstream of shock boundary layer interactions 
with bleed. 

Hamed and Shsng [1] reviewed the existing 
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experimental data base for shock wavc/boundary layer 
interactions in supersonic inlets and other related 
configurations. According to this survey, there is enough 
experimental evidesoe P-14] to indicate that local bleed 
-- can control flow separation in shock wavc/boundary layer 
interactions. There are disagreements [1] however among 
the different experimental results regarding the effects of 
bleed location relative to the shock P-12]. Strike and 
Rippy P] measured the surface pressure in the interaction 
zone of an oblique shock wave impinging a turbulent 
boundary layer over a flat plate with suction. They 
determined that less suction is required to control 
sepa r ation, when applied upstream of the shock. Sccbaugh 
and Childs [10] investigated experimentally the 
axisymmctric flow in the interaction region of the boundary 
layer inside a duct. Contrary to the conclusions of Strike 
and Rippy P], suction within the interaction region was 
found to be more effective in suppressing the effects of 
separation. 

Remlinger et aL [IS] and Hamed ct aL [16-20] 
modeled bleed by resolving the flow through ibe individual 
bleed ports. The advantage of this approach ever that 
which si mul ates the bole roughness effects [13, 14] is that 
the bleed mass flux distribution, the bleed drag, as well ss 
the boundary-layer characteristics downstream can be 
determined from the numerical solution. Comparative 
information regarding the effect of bleed hole sizes, slant 
angle and bleed location relative to the shock can be 
obtained using this spproachto help resolve the o oo trov eny 
among the different experimental data concerning the effect 
of these parameters. Hamed et aL [17] investigated the 
effects of slot width, depth [18-20] and slot stoat angle on 
the flow field in the inter a ction zone of an oblique 
• shock/turbulent boundary layer interaction. The raeuks 
obtained in these investigations for bleed through slob 
located across the shock impingement point indicate that 
turbulent boundary layer separation eould be eliminated by 
removing 10 to 16 percent of the incoming boundary layer 
mass flow through slanted slots. The slots with low stoat 
angle have higher discharge ooefficieots and start to bleed 
at higher plenum pressures. In addition they produce 
redeveloping boundary layer downstream with higher 
friction co, ffletent. The purpose of this study to to 
investigate the effect of bleed location relative to the shock 
on the flow field in the interaction region aod downstream. 
For this purpose numerical simulations were conducted for 
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both normal and 20* slanted slots at three locations, 
tqiatream, across and downstream of the shock 
impingement point. The computed results for the different 
configurations are presented to demonstrate their effect on 
the flow field in the interaction zone and downstream of the 
interaction as well as the bleed flow inside the slot. 

Commits tkmal f*-th 

Numerical flow simulations were conducted for 
bleed applied through normal and slots in the 

interaction region of an externally generated incident 
oblique shock and a turbulent boundary layer. The flow 
was considered to be turbulent throughout the calculation 
domain; no attempt was made to model transition and/or 
relaminarization. As in our previous studies [16-20], die 
solution domain used in the two dimensional flow 
simulations extends upstream of the fist plate leading ed g e , 
downstream of th: interaction region and inside the bleed 
slot, where the specified plenum pressure at die bottom 
controlled the amount of bleed mass flow. 

The upper solution domain boundary was set at 
0.13 ft. from the flat plate with the lower solution domain 
at the bottom of the bleed slot. In the streamwise 
direction, the domain extended along the plate length of 
1.13 ft. and the distance of 0.07 ft., upstream of the plate _ 

leading edge to the station where the uniform u p stream 
conditions are applied. The incident oblique shock crosses 
the upper boundary, and all other wave systems including 
the reflected, and any separation and roattaohmsat shocks, 
cross the outflow boundary. Uniform freeatream flow 
conditions are applied at die inflow boundary and all 
va ri a bles are ext rapol ated at the outflow boundary. No alp 
a di a b a tic flow conditions are applied at the plate and bleed 
walla. The static pressure is specified and first order 
extrapolation is applied for the rest of the flow variables st 
the bottom of the bleed slot The location of the 
•hock is fixed at the upper boundary with free stream end 
post shock conditions specified upstream and downstream. 

The PARC code [21] with a two-equation kt 
turbulence model was used in this bleed study. The 
turbulence model is based on Chien’s low Reynolds number 
k-« model [22] with Nichols [23] modifications to add 
compressibility effects. This code has been used in several 
supersonic and hypersonic inlet flow computations [7, 24]. 

The computations were performed using a 308x68 
grid over the plate surface and an 89x68 grid ioxtto the 
bleed slot Variable grid spacing was used in both x and 
y directions for grid clustering around the bleed walk and 
at the plate curfoee, with y*. - x_* « 0.2579410 4 ft. 
oorroaponding to y* - 2.0 at x ■ 0.9 ft. The 
computations for a typical bleed ease required 9000 local 
time stops at 0.2 CFL number to reach Heady state sohttfon 
baaed on six orders of magnitude reductions in the 
averaged root mean aquare error in the flux at three 
different locations relative to the shock. 


The computations were conducted at the incoming 
flow conditions of M« 2.96, Re « 1.2 x 10 7 /ft, and an 
impinging oblique shock angle of 29.84* (wedge angle of 
7 .93*) corresponding to the tost conditions of reference [29] 
for the separated flew case with no bleed. The flow 
characteristics inside the bleed slot, and throughout the 
interaction were determined for bleed through a normal and 
20* slanted slot. The slot width for the normal bleed case 
was 0.01069 ft which is equal to 0.8085 tones the 
boundary layer thickness upstream of the interaction. To 
maintain comparable bleed mass flows among the normal 
and slant: d bleed cases, the opening at the plate surface 
was maintained at 0.0106 ft leading to a slot width which 
is 0.2766 times the boundary layer upstream fin the 20* 
slanted slots. The effect of bleed location was investigated 
by studying the flow field with bleed applied upstream, 
•cross and downstream of the incident abode interaction 
with the plate surface. The slot opening center, 3^, was 
located at 0.99468, 1.0 and 1.00532 ft. from the plate 
leading edge for the three cases capeetivcty. This 
corresponds to the incident shock interaction point with the 
plate was situated at the downstream corner, the center and 
the upstream comer respectively for the three bleed 
locations. 

Figure l presents the Mach number contours for 
the six different slot geometries at choked conditions. The 
wave system across the Meed opening consists of an 
expansion fan at the upstream slot comer then en oblique 
■bock at the downstream comer of the darted bleed slot or 
a bow shock that originates inside the normal slot. The 
"circulation flow region inside the slot at upstream wall 
occupies a large portion of the slot width in the nesmal 
bleed cases, but is much smaller in extent ; "*“ u the 20* 
a l a n t ed slot. The two mechanisms controlling the 
formation of recirculation region inside the slot at the 
upst r ea m wall are insufficient flow turning into the normal 
bleed (lot, and the impinging oblique abode wave which 
originate* at the downstream edge of the bleed slot in the 
20* slanted bleed [20]. The corresponding Mach number 
and flow angle distribution across the skit opening a re 
presented in Figures 2 and 3. The static pressure and 
friction coefficient distribution over the plate surface are 
* presented in Figures 4 and 5. In all these figures, the 
results for toe three different bleed locations ween 
such that they were all aligned around the center of the 
bleed opening, in each ease. According to these 
figures, the Mach number kinkklly higher in the upstream 
bleed case but approaches the same value towards tb* slot's 
downstream oomer as the other two bleed locations. The 
flow angle k equal to the slot angle k the 20* slanted Masd 
but k less than 40* in the oaae of normal bleed. In the 

normal bleed slot eases, the flow angles am highest for the 
downstream bleed and lowest for the upstream bleed. Urn 
initial (harp flow overturning in the normal bleed ease 
(Fig. 3) k caused by a small separation bubble on the pkto 
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surface upstream of the (lot which it confirmed by the 
negative friction in Fig. 5. According to Fig. 5 which 
preaenta the skin friction diatribudon over the plate surface, 
a amall icparated flow region ia predicted at upstream of 
the (lot and an even amalter one at downatream comer in 
all normal bleed caact. On the other hand, no flow 
reparation wu predicted over the plate turface with the 20T 
Manted bleed (lot upatream and acrota of the abode. In the 
20* alanted bleed downatream of the (hock, the flow 
aeparatea over the plate aurface upatream of the (lot 
opening. According to Figs. 4 and 5, the maaa removal 
effect! extend both upatream and downatream of the bleed 
alot and ia reflected in both the aurface preaaure and 
friction coefficient. The preaaure rite downatream of the 
bleed (lot ia faster in the normal bleed than the alanted 
bleed for 1 three locations. The change in die friction 
coefficient over the plate aurface downatream of die 
interaction ia characteriaticaUy that of boundary kyer 
redeveloping in the eaae of bleed through the aknlrid alota 
[20]. However the akin friction distribution downatream of 
the interaction ia very different in the eaae of normal bleed 
since the stagnation point is slightly inside the (lot’s 
downstream wall [20]. 

Figure 6 presents the variation ia die discharge 
coefficient with the plenum pressure, P k , normalized by the 
local pressure, P*. For die bleed locations upstream, 
across and downstream, P*. represents the pressure before 
the incident shock (P*. <* 1.0 P^), behind the incident 
shock (Ph, ” 1.7752 P* ), and behind the reflected shock 
(P ta * 2.96 Pw) respectively. The same data ia 
represented in Fig. 7 in the form of bleed mass flow as a 
perccntageof the boundary layer mass flow upstream of the 
interaction at x ■ 0.9 ft According to these figures, the 
highest discharge coefficient is obtained with the upstream 
bleed through 20* alot and the lowest with downstream 
through normal afc*. The alanted bleed discharge 
coefficient is very sensitive to the bleed location where 
higher discharge coefficients are obtained with upatream 
bleed, while .«e downatream bleed resulted in c hok ed 
discharge coefficients lower than the across the s hoc k 
bleed. Tbe difference between the discharge coefficients 
for the upstream and across the shock normal bleed is very 
small. The highest bleed maaa flow U obtained with 
normal bleed across the shock and with 20* alanted bleed 
downstream of tbe shock. The lowest choked bleed mass 
flow is obtained with the normal bleed u p st r ea m of the 
ahook. In general bleed mass flow changes were gradual 
and choking oeeurxed at lower plenum pressures in the case 
of normal bleed. This is consistent with tbe experimental 
hole blood data presented by Sybcrg and Koncack [2]. 

Tbe effect of bleed mass flow on tbe boundary 
kyer characteristics downstream of tbe int era ctions at 
X»1.13 ft k pr esente d in Figures I through 11 for all (fat 
configurations. According to these figures, all slot 
configurations produce nearly the same dk p kocmcut 
thickness, 3* ■ 0.0024 ft. and momentum thi c kncee , t ■ 
0.0007 ft. at choking, but the bleed through the normal slot 


across the shock produces lower values (5* » 0.0020 ft. 
and 9 ■ 0.0006 ft.) at a tower bleed mass flow rate (8* of 
the boundary kyer). In general the transformed form 
factor downstream of the interaction did not change 
aignificantly over tbe bleed mass flow range. It had the 
smallest value (H* ■ 1.36) for the normal alot at 15% 
bleed mass flow. This is consistent with the data presented 
by Fukuda ct al. [3] for the scoop in the case of the 
forward ccnteibody interaction where the upstream edge 
Mach number wu 2.09. However the friction coefficient 
is low for the normal bleed at these coodkioas. Tbe small 
effect of the bleed configuration can be attributed to the 
small differences in the bleed slot locations (a maximum of 
0.0106 ft. or 0.8085 daws the boundary kyer thioknua 
upstream). Additional results were obtained for wider 
normal slots, whose centers were located at X“ 0.97872 
ft., X»1.Q ft and X»1.04256 ft. from the leading edge. 
The resulting boundary kyer characteristics for these wide 
slots (D/5 «• 3. 234) are presented in Figures 12 and 13. 
According to the results in these figures, across the shock 
bleed produces larger reductions in both disp l a cem ent and 
momentum thickness than both upstream and downstream 
bleeds. AU these configuration produce the same effect at 
bleed mass flows equal to 28% of the inooming boundary 
kyer. However these large bleed mass flow with the wider 
normal slots do not improve tbe downstream boundary 
kyer characteristics compered to the smaller bleed mass 
flows through the narrow alota. 

Conclusions 

A numerical investigation wu conducted to study 
the effect of 20* skated and normal bleed alot location 
relative to the shock on tbe flow field in an oblique shook 
wave/turbulent boundary kyer interaction. The results 
indicate that the slanted bleed alot discharge coefficient is 
more sensitive to bleed location. The highest discharge 
coefficient k obtained with the skated bleed alot located 
upstream of the incident shock followed by the alot located 
across the shock. The normal bleed discharge coefficient 
k leu sensitive to the bleed location but k lowest ia dm 
oaae of downstream bleed. Higher bleed mau flows are 
obtained boweverwith tbe downstream blood location. The 
boundary kyer momentum and di apkeem e n t thickne s s 
downatream of the interaction approached the same value 
for all bleed configurations at choked co ndi tio n s, but 
• otherwise varied with the bleed location in tbe caw of 
normal unehoked bleed. The amaflaat reduction in the 4* 
and $ were obtained with normal upstream bleed and dm 
krgeat reduction with dm across dm shock un e hoked 
normal bleed. However tbe friction coefficient wu much 
higher for alanted bleed. The boundary kyer 
characteristics Were more sensitive to bleed location for 
larger bleed mau flows through wide normal alota. In 
general aerou the shock bleed* resulted in larger 
reductions in both f and # but all bleed location* war* 
moat effective and produced similar reductions when 26 
percent of the inooming boundary kyer wu bled through 
the normal alota. 
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Abstract 

Numerical simulations were conducted to 
investigate the performance characteristics of bleed 
through normal slots and its effect on the turbulent 
boundary layer development under zero and strong adverse 
pressure gradient caused by incident oblique shock. The 
solution to the compressible Navier-Stokes and k-E 
equations was obtained in a domain that includes the 
regions inside the bleed slot and plenum in addition to the 
external flow. The computational results demonstrate the . 
interactions between die plenum, and bleed flow and die 
effect of incident shock on the boundary layer development 
downstream. The computed results agree with the 
experimentally measured pitot and static pressure 
distribution inside the slot. The bleed mass flow without 
incident shock was underpredicted over the range of 
plenum pressures. The computations predicted die 
measured increase in bleed mass flow with incident shock. 

Introduction 

Bleed is used in supersonic inlets to control the 
effects of flow separation associated with shock boundary 
layer interactions on ramps, cowls and side walls. In order 
to predict the effect of bleed on the boundary layer 
development and internal shock structure, few investigators 
[1-3] included ramp chambers, throat plenum and exit 
louvers in their supersonic inlet flow simulations. 
However, because of the enormous increase in 
computational time and grid generation complexity, most 
investigators simulated the global effects of bleed by 
changing the boundary conditions [4-7] and/or die 
turbulence models [8,9] in the bleed regions. The imposed 
mass flux in the bleed area was based on experimental 
measurements in some investigations [7] and on empirical 
correlation of the bleed discharge coefficient in others [4- 
6]. Even when the turbulence model was modified for 
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mass removal at the porous wall [8] the effect of bleed on 
boundary layer development was not adequately predicted. 

Reproducing the experimentally measured boundary layer 
profiles over the porous wall even required the 
implementation of blowing rather than suction in the 
turbulence model [8]. Paynter et al. [9] used a rough wall 
algebraic turbulence model to simulate the increase in the 
growth rate of boundary layer induced by bleed, by 
changing the roughness parameter in the numerical 
solutions to the Navier-Stokes equations until the boundary 
layer profiles matched the experimental measurements [10, 
1 1]. They determined that roughness was a strong function 
of the upstream boundary layer mass flux removed, in the 
bleed rates between 3% and 15%. However Paynter et al. 
pointed out that tire experimental set-up whose data they 
used to develop their model, may have produced non 
uniformity in the mean velocity distribution and 
nonequilibirum boundary' layer in the bleed region. In 
general, bleed models are restricted by the particular flow 
conditions and any inadequacies in the experimental data 
used in their model development. According to Hamed 
and Shang [12], the conclusions regarding the effects of 
bleed hole size and bleed location relative to the shock 
were not consistent among several experimental studies. 

Hamed, et al. [13-18] and Remlinger, Shih and 
Chyll [19-22] followed an alternative approach to 
investigating bleed effects. They conducted numerical 
simulations in which the viscous flow field was resolved 
inside the individual bleed holes and slots. This approach 
helps in understanding the important local phenomena that 
control the flow in the bleed regions, and how they are 
affected by the bleed configuration, and the external flow 
conditions. These type of investigations can and have been 
used in simple flow configurations to study the relative 
effects of bleed/slot angle, size and location. They 
revealed the presence of a separation bubble near the 
hole/slot entrance and a "bleed shock" that initiates inside 
the hole/slot and attaches to the downstream comer under 
certain conditions. The size of the separation bubble which 
controls the bleed discharge coefficient is strongly affected 
by bleed hole/slot slant angle and less by the free stream 
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conditions and plenum pressure. Tie boundary layer 
development downstream was strongly influenced by the 
bleed angle and the bleed shock strength. 

Hie purpose of present numerical study is to 
investigate ihe effects of external flow and plenum 
conditions on’ bleed performance. In the investigated 
configuration bleed is applied to a flat plate turbulent 
boundary layer through a normal slot. The numerical 
solution to die viscous flow inside the slot and plenum and . 
over the plate surface were obtained for supersonic flow 
with and without impinging shock at the slot's upstream 
comer. The flow computations were conducted over a 
range of plenum pressures up to choked bleed conditions. 
The results, which are compared with experimental data, 
show the flow characteristics inside the bleed slot and die . 
effect of bleed on the boundary layer development 
downstream with and without incident shock. 

Flow Configuration a nd Computational Details 

In the investigated flow configuration shown 
schematically in Fig. 1, bleed is applied through a normal 
slot to the turbulent boundary layer on a flat plate. The 
bleed slot whose width is one centimeter and its depth 2.S4 
cm is connected to a 41.4 cm x 63.S cm plenum. Outside 
of the slot, the solution domain extended 12.2 cm above ., 
and 40.9 cm along the flat plate surface. The boundaries 
AB and CD Were located at 28.2 cm upstream and 12.7 cm 
downstream of the slot's upstream comer to match the 
locations of experimental velocity profile measurements in. 
reference [6]. 

The flow computations were conducted using the 
PARC code [23] with the compressible k-e turbulence 
model of Chien [24] modified by Nichols [25] for 
compressibility effects. The change in bleed mass flow 
rate was achieved through changing die static pressure at 
the plenum's outflow boundary. Referring to Fig 1, the 
upstream boundary conditions consisted of free stream 
Mach number, stagnation pressure and stagnation 
temperature and a computed boundary layer velocity 
profile for flat plate computations which was matched to 
the experimentally measured profile's displacement and 
momentum thickness. The pre-shock and post-shock 
conditions were specified along the upper boundary to 
locate the inviscid incident shock at the upstream comer of 
die slot and flat plate surface. Flow variables were 
extrapolated at the downstream boundary and at the 
plenum's outflow boundary where the static pressure was 
specified to control the bleed mass flow. The initial 
conditions were specified from the solution without bleed 
for a given free stream Mach number. Subsequently, die 
computed flow field with the bleed slot and plenum was 
used to initialize the solutions at higher bleed mass flow 
rates (lower plenum pressure settings). 


The computational grid shown in Fig. 2 includes a 
308x68 grid over the flat plate surface, 89x68 grid inside 
the bleed slot and 187x148 grid inside die plenum 
chamber. Variable grid spacing was used in both x and y 
directions for grid clustering around the bleed walls, 
plenum chamber walls and at the plate surface, with 

x . = 0.9446 xlO -3 cm and y . = 0.4755 x 10~ 3 cm 
min * mm 

corresponding to y + = 0.959 at the inflow turbulent 

boundary layer at 2.46 free stream Mach number. 

Results and Discussion 

Typical results of the computed flow inside the 
slot near choked conditions are presented in Fig. 3 and. 
compared with the experimental results of Davis et al. [27] 
reproduced in Figures 4 for incident oblique shocks whose 
inviscid impact point coincides with die slot's upstream 
comer at x=0. The figures present die pitot and static 
pressure contours for a shock generation angle a - 8°, at 
2.46 free stream Mach number. The computations predict 
very well both the shape and magnitude of die pressure 
contours inside the slot, the location of flow separation, 
and the size of the separation bubbles on die slot walls. 
The computed Mach number and velocity vectors, which 
are shown in Fig. 5, indicate flow angles of 33° from the 
plate surface, and high Mach numbers of 1 .8 in a triangular 
region at the center of the slot opening. Disagreements can 
be observed between die experimental and computational 
results in this region. Since the pitot probe was always 
parallel to the slot walls in the experiment, one should 
expect discrepancies in these regions where the probe is 
oriented at 67° angle to the Mach 1.8 flow. The sharp 
gradients in the pressure across "bleed shock" were also not 
resolved by the measurements. The experiment indicates 
reattachment of the flow on the slot's downstream wall near 
the exit, which was not predicted by the computations. 

Effect of Incident Shock 

The effect of incident shock at the slot's upstream 
comer is demonstrated by comparing die computational 
results with and without the shock. Figures 6 and 7 present 
the computed Mach number and pressure contours inside 
the slot and in the neighboring regions over die plate 
surface and inside the plenum. The Mach number contours 
inside the slot indicate that when there is no incident shock, 
the separated flow region along the slot's upstream wall is 
much larger while die separated flow region over the slot's 
downstream wall near the exit is much smaller. In 
addition, the Mach number contours indicate less flow 
turning into the slot when there is no incident Shock. This 
coincides with a weaker expansion fan at the slot's 
upstream comer as can be seen from the pressure contours. 
Also without incident shock, the "bleed shock" which 
originates inside the slot is quickly weakened outside 
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through interactions with a second expansion fan at the 
slot's downstream corner. 

The maximum computed bleed mass flow with 
and without incident shock was S.1% and 2.35% of the 
incoming boundary layer for the slot geometry (D/6 « 
0.38). The pressure and Mach number contours indicate a 
separation bubble on the plate surface upstream of the slot 
in the case of incident shock. This is confirmed by the skin 
friction coefficient distribution of Fig. 8. One can see from 
Fig. 8 that there is a large difference in die predicted 
friction coefficients with and without incident shock within 
two slot widths downstream of bleed. In order to 
understand the cause of this difference, the mass flux, flow 
angle, and Mach number distribution across the slot 
opening are presented in Figs. 9-11. In the case of no 
incident shock, the mass flows out of the slot in the 
subsonic region behind the "bleed shock" and the 
downstream slot wall. This flow turns around die slot 
comer and forms die boundary layer over the plate surface. 
In the case of incident shock the subsonic flow goes into 
the slot opening in the region behind the "bleed shock", 
and the flow over the plate surface downstream originates 
outside the slot. 

According to the velocity profiles of Fig. 12, the 
effect of incident shock is seen to extend over two to three 
slot widths upstream and downstream of the slot, where the 
distorted velocity profiles are less full near the plate. The 
distortion upstream is associated with the shock induced 
pressure gradient and flow separation. The shape of the 
profiles downstream is affected both by the local flow 
conditions above the slot and by the strength of die "bleed 
shock" at the slot opening. Since the incident shock 
reduces the flow velocity, the profiles downstream of the 
slot are less full near the wall. The "bleed shock", has 
smaller effect on the downstream velocity profiles hi the 
case with no incident shock because it is quickly weakened 
by the expansion fan at the slot's downstream comer.. 

Plenum Interactions 

A considerable part of die computational effort 
was consumed in modeling the flow field inside the 
plenum. The increased computational effort was not only 
associated with the number of grid points inside the 
plenum (50.6% of the total) but also with the slower 
convergence of the flow at low subsonic plenum speeds. 
In order to determine the nature of die plenum interactions 
with the bleed and external flow fields, the bleed mass flux 
was monitored at the slot opening (bleed inflow) and exit 
(bleed outflow) during the computations. Fluctuations 
persisted in the bleed mass at both inflow and outflow 
boundaries when die computations were conducted using 
local time stepping. The amplitude of the bleed mass flow 
fluctuations was higher at higher plenum pressures and 


diminished near choking. Subsequendy one time accurate 
flow simulation was performed with o*8° incident shock 
at choked bleed conditions. Ihe computational results 
indicated no fluctuations in the bleed mass flow at the slot's 
inflow boundary in this case and much smaller fluctuations 
in the bleed mass flow at the outflow boundary compared 
to die local time stepping results. Additional numerical 
computations with local time stepping were then performed 
without the plenum. The bleed mass flow in these 
computations did not fluctuate and coincided with die time 
accurate bleed mass flow at the slot’s inflow boundary. 

The discharge coefficient, Q, and the flow 
coefficient, F, are used to present die bleed mass flow data 
in nondimensional form. They are defined as follows: 

1 T+1 

Q = m b / A b PA-l—Y (— ) 2 <’-» 

* 4 4 r RTj y+r 

When the stagnation free stream value P is 

l iaf 

used for the reference pressure, P p die discharge 
coefficient, Q, represents the ratio between the bleed mass 
flow and the ideal mass flow that could pass into the slot at 
sonic conditions if the flow was to expand isentropically to 
fill the bleed area A b . On die other hand, the flow 
coefficient represents the ratio between die bleed mass 
flow and the mass flow that would pass through the bleed 
area at the flow conditions p c U e outside the boundary 
layer 

F = m b / 

The predicted variation in die discharge 
coefficient with plenum pressure is presented in Fig. 13 
with the experimental results of references [26] and [27]. 
The figure shows the computed bleed mass flow without 
plenum, as well as the range of bleed mass flow 
fluctuations in the computations with plenum using local 
time stepping. The computations without plenum agree 
with die experimental data in die case of incident shock, 
but under-predict die bleed mass flow without shock over 
the range of plenum pressures. The much higher discharge 
coefficient in the case of incident oblique shock can be 
caused by the increased static pressure above die slot 
and/or the velocity distortion upstream. In order to 
separate the two effects, a bleed discharge coefficient 
based on the local static pressure, P^, is presented in Fig. 
14. The value of the local static pressure was taken as P^f 
for no incident shock and Pj, the inviscid static pressure 
behind die incident shock (Pj/Pinf “ 1.7). This figure 
indicates that the difference in the computed discharge 
coefficient with and without incident shock becomes 
smaller based on this normalization. The remaining 


difference can be attributed to the incoming velocity 
profile distortion by the incident shock. 

Conclusions 

Numerical simulations were conducted for 
supersonic flow at Marh 2.46, with bleed applied through a 
normal slot with and without incident shock at die slot's 
upstream comer. The computed results agree with the 
experimental measurements inside the slot in the case of 
incident shock. The computations indicate higher flow 
angles at the slot opening and a narrower separation bubble 
on the slot’s upstream wall in the case of incident shock 
which leads to higher discharge coefficient. The "bleed 
shock" which originates inside the slot and extends as -an 
oblique shock outside is quickly weakened in die case with 
no incident shock by the interactions with an expansion fan 
at the slot’s downstream comer. 

The plenum caused fluctuations in the bleed mass 
flux at the slot exit but time accurate computations at 
choked conditions indicated that these distuibances did not 
affect the flow and mass flux at the slot opening. The 
fluctuations were associated with the lateral movement of 
the bleed flow between die two separation bubbles on the 
slot walls. The computed bleed mass flow without plenum 
coincided with die predicted steady value at die slot inlet in 
the time accurate numerical simulations with plenum. The 
experimental data for bleed mass flow agreed with the 
computations in the case of incident shock but were higher 
with no incident shock. 
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Fig. 1. Schematic of Com putational Domain. 


Fig. 2. Computational Grid Near Bleed Slot, 
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Fig. 5. Velocity Field Inside Slot 




(a) With Incident Shock 

Fig. 6. Comparison of the Mach Number Contours With 
and Without Incident Shock 
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(a) With Incident Shock 


(b) Without Incident Shock 


Fig. 7. Comparison of the Pressure Contours With and 
Without Incident Shock 
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Fig. 8. Skin Friction on die Plate Surface. 


Fig. 9. Bleed Mass Flux Across Slot Opening 
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Fig. 10. Flow Angle Across Slot Opening. 


Fig. 11. Mach Number Across Slot Opening. 
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(a) Upstream of Bleed Slot (b) 

Fig. 12. Boundary Layer Velocity Proxies 


(b) Downstream of Bleed Slot 
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□aeon Experimental Data 
***** Computation w/o Plenum 
Computation w/ Plenum 





